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We present a detailed study of a monolayer film of 3,4,9, 10-perylene-tetracarboxylic acid dian- 
hydride (PTCDA) on Ag(lll) (the PTCDA/Ag(lll) system). The study is done within density 
functional theory with the use of the periodic slab model. The slab is chosen to contain a PTCDA 
monolayer film on a silver thin film of different thicknesses (6, 9, and 12 layers) with the (111) 
orientation. We show that one of two surface states of the pure Ag(lll) films transforms into an 
unoccupied interface state due to the adsorbate-substrate interaction. The relation of the resulting 
state to the unoccupied state that has been experimentally observed in the PTCDA/Ag(lll) system 
by scanning tunneling and two photon photoemission spectroscopy is discussed. 

PACS numbers: 73.20.-r, 73.20.At 



I. INTRODUCTION 

Theoretical and experimental investigations of elec- 
tronic properties of organic molecular structures ad- 
sorbed on metal surfaces are of great importance for the 
solution of both fundamental and practical problemsii^. 
Among organic molecules, the particular research inter- 
est is expressed in 3,4,9, 10-perylene-tetracarboxylic acid 
dianhydride (PTCDA) molecules, since the latters are 
able to form strictly ordered layers on the surfaces of 
some metals. ^^"^ When adsorbed on Ag(lll), the PTCDA 
monolayer film causes the appearance of a delocalized 
two-dimensional band state that in Ref. ^ has been re- 
vealed at ~ 0.7 eV above the Fermi level within scanning 
tunneling spectroscopy (STS) observations. This unoc- 
cupied state is described by a parabolic dependence on 
the two-dimensional electron wave vector k with the ef- 
fective mass m* = 0.47me, where rrie is the free electron 
mass. Aiming at understanding the mechanism of elec- 
tron delocalization in molecular structures, in spite of the 
similarity of m* with that of the Shockley surface state 
(SS) of the pure Ag(lll) surface, the authors of the cited 
work has attributed the origin of such a dispersing state 
to the substrate- mediated interaction of LUMO+1 or- 
bitals of ordered arrays of the PTCDA molecules. The 
main experimental evidence of it is that the spatial dis- 
tribution of the state is concentrated on the molecules 
and resembles that of the LUMO+1 level of PTCDA. At 
that the anhydride groups of PTCDA are involved in the 
wave function of the state. ^ 

In Ref. 6, the unoccupied state of the PTCDA/Ag(lll) 
system has been investigated by means of two-photon 
photoemission spectroscopy (2PPES). The authors of the 
cited work have presented experimental evidence that the 
state in question (here with the energy of 0.6 eV above 



the Fermi level at k = and with a parabolic disper- 
sion with m* = 0.39 ± 0.03me) comes from the SS of 
Ag(lll). The argumentation in favor of such an origin 
of the state is that their angle-resolved 2PPES observa- 
tions reveal the effective electron mass identical to the 
SS of the pure Ag(lll) surface and the location of the 
state at the metal-organic interface with large bulk pen- 
etration comparable to that of the Shockley state. Thus, 
there are at least two diametrically opposite points of 
view on the origin of the dispersing unoccupied state 
in the PTCDA/ Ag(l 11) system. Since such states play 
a crucial role in the dynamics of excited electrons and 
charge transfer through the interface^^^, the question of 
the origin and properties of the unoccupied state in the 
PTCDA/ Ag(l 11) system is important to be cleared up. 



In this paper, within the periodic slab model we study 
the electronic structure of the PTCDA/ Ag(l 11) system, 
using ab initio density functional theory calculations. 
The slab is chosen to contain the PTCDA monolayer film 
on Ag(lll) films of different thickness (6, 9, and 12 lay- 
ers). In the monolayer, the PTCDA molecules are packed 
in a herringbone structure. On the one hand, we show 
that the experimentally observed unoccupied state arises 
from the Shockley surface state, has the bulk penetra- 
tion similar to that of the latter state, and is localized 
at the interface. On the other hand, we demonstrate 
that in the plane of the PTCDA monolayer film the spa- 
tial distribution of this interface state is concentrated on 
the molecules with quite big amplitude on the perylene 
backbone and anhydride groups and relatively small - 
at the perylene edges of the PTCDA molecule. Thus, 
we corroborate the interpretation proposed in Ref. 6 and 



explain the STS observations of Ref. 
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FIG. 1: (Color online) Structure of the PTCDA monolayer 
on the Ag(lll) surface. 



II. CALCULATION DETAILS 

We use the periodic slab model to study the metal- 
organic interface. The electronic structure calcula- 
tions were performed by using the OPENMX code^^ 
which are based on the density functional theory (DFT) 
and a linear combination of localized pseudoatomic 
orbitals method. We exploit Troullier-Martins type 
pseudopotentials^^^ with a partial core correction to re- 
place the deep core potentials by norm-conserving soft 
potentials. Upper core electrons for the Ag pseudopoten- 
tial are included in addition to valence electrons in order 
to take into account the contribution of semicore states to 
the electronic structure. As a basis set, we take double- 
valence plus single-polarization orbitals. For hydrogen, 
carbon and oxygen atoms constituting PTCDA, varia- 
tionally optimized basis orbitals for biological molecules 
are used with cut-off radii of 5.0 for C and 4.5 a.u. for 
O and H atoms. For Ag atoms cut-off radius is 7.5 a.u. 
The real-space grid techniques are used with the energy 
cutoff of 190 Ry in numerical integrations and solution 
of the Poisson equation with the help of the fast Fourier 
transformation. The total-energy convergence was better 
than 0.027 meV. 

To describe the PTCDA/Ag(lll) system, we employ 
a slab containing a silver film with the (111) orienta- 
tion and Ni atomic layers (A/'^=6, 9, 12) together with a 
PTCDA monolayer film, which is "applied" on one side of 
the silver film only, and a vacuum region corresponding 
to 6 silver interlayer spacings. The distance between the 
silver surface layer and the PTCDA monolayer is taken to 
be equal to experimental value of 2.86 A (see Refs. 13 and 
fl3 ). A close value can be obtained in DFT calculations^^ 
with exchange-correlation potential in the local density 
approximation (LDA), while in the generalized gradient 
approximation the distance becomes very large (~ 3.5 A 
in Ref. tl3i and .16.) and the adsorption energy is under- 



estimated. This fact causes the use of the LDA in our 
DFT calculations. 

The PTCDA monolayer on Ag(lll) has an ordered pe- 
riodic structure that consists of two types (denoted as A 
and B in Fig. [1]) of molecules misoriented by an angle of 
77° with respect to each other (a herringbone pattern). 
The A-type molecule is aligned parallel to the row of sil- 
ver atoms, while the B-type molecule is misaligned by an 
angle of 17° (see Refs. ITtI and HHI) . The superstructure of 
the interface is given by a nearly rectangle (a deviation 
of 1°) unit cell of 18.96 x 12.61 A^ The resulting peri- 
odic supercell of the slab model contains 33 x Ni silver 
atoms (33 - the number of atoms per layer) and 76 atoms 
belonging to two PTCDA molecules. The Brillouin zone 
of the supercell is sampled with a (3 x 3 x 1) mesh of k 
points. The z-axis is directed normal to the surface of 
the films. 



III. RESULTS AND DISCUSSION 

We start with an band-structure analysis of clean 
Ag(lll) films with the supercell that contains 33 x Ni 
silver atoms in order to identify surface states (SS) in 
the absence of the energy gap at the F point of the corre- 
sponding small surface Brillouin zone and to make sure 
that characteristics of these states are the same as in 
the case of the supercell with 1 x Ni silver atoms, which 
gives the familiar surface band structure for the Ag(lll) 
surface. Note that in the case of a film there are two 
surface states with different energies caused by bonding- 
antibonding splitting due to interaction through the film. 
We distinguish these two states as SSI and SS2 with en- 
ergies Essi > Ess2 at the f point. A comparison has 
shown that the surface-state energies obtained with the 
33 X A^^ supercell (see Table |T]) coincide with those found 
with the 1 X Ni supercell (not presented) within ~ 10 
meV, while the bonding- antibonding splitting {AEdean) 
— within ~ 2 meV. As regards the effective masses, their 
values practically do not vary with the number of atoms 
per layer. 

The mentioned characteristics of the surface states de- 
pend on thickness of the film (see Table [I]). As was ex- 
pected, AEciean dccrcascs rather rapidly with increasing 
Ni and should vanish at the Ni ^ oo limit. The en- 
ergy Essi tends to reduce, while Ess2 seeks to increase. 
At the same time the effective mass m^g-^ varies weakly 
with the thickness as compared to ^552- Note that as 
Ni increases the average energy ^^^Zean becomes closer to 
the surface-state energy E^^^ = — 63±1 meV experimen- 
tally observed^^ for the Ag(lll) surface and already for 
the twelve-layer film ^^Zean differs from Eg^g within ~ 16 
meV. 

Now we consider the PTCDA monolayer on the 
Ag(lll) film of different thicknesses. Fig. [2] shows the 
obtained band structure for two cases: nine- and twelve- 
layer silver films. Also, the figure reflects our results on 
analysis of the square modulus of wave functions of found 



0.2 0.15 0.1 0.05 0.0 0.05 0.1 0.15 



0.2 0.15 0.1 0.05 0.0 0.05 0.1 0.15 



FIG. 2: (Color online) Band structure of the PTCDA monolayer on the nine- (a) and twelve-layer (b) Ag(lll) film in the 
X2 ^ F ^ Xi direction of the surface Brillouin zone. The thick gray lines show the parabolic approximation of the dispersion 
of the interface (upper line) and surface (lower line) states. Points indicate wave vectors k, at which the spatial distribution of 
the interface-state wave function is analyzed (see Table [III and Fig. [S). 



TABLE I: Energies Essi and Ess2 of the surface states at the 
r point and corresponding effective masses mgsi and mgs2 
for clean Ag(lll) films with different number of layers (Ni). 

AEclean = EsSl - EsS2 and ^,7,,^ = EsS2 + /^Eclean/2. 

In the case of the PTCDA monolayer on the Ag(lll) films, 
Eis and Ess are the F-point energies of the interface and 
surface state with effective masses mj^ and mj^, respectively. 

= Eis - Ess and Eav = Ess + Ah the energies 

are measured electron volts from the Fermi energy. 
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states (see the thick lines). When averaged over spatial 
coordinates xy at different k along f ^ X2 and f ^ Xi 
directions, the square modulus can be considered as a 
function of which gives important information about 
charge-density distribution [denote as p(z,k)] along the 
direction normal to the surface. We have thus revealed 
the partly occupied surface state (the lower thick line in 
Fig. [2]) with the wave function localized in the region of 
the "clean" surface of the Ag(lll) film. At the F point, 



the energy (denote as Ess) of such a state depends on the 
film thickness and, as is seen from Table in tends to E"^^^^^ 
with increasing Ni. This surface state is originated from 
the state marked as SS2 in the case of the clean Ag(lll) 
films. Note that in spite of the difference Ess — Ess2, 
which is caused by the presence of the PTCDA monolayer 
and decreases with an increase of Ni (see Table [J), the 
effective masses and TTi ^ are the same (except for 
the six-layer silver film). At that, for the twelve-layer sil- 
ver film they become very close to the experimental value 
of 0.40me for the surface-state effective mass found for 
the Ag(lll) surface from photoemission data in Ref. [isl. 

On analyzing p(2:,k), also we can unambiguously say 
that in the spectrum above the Fermi level there is an in- 
terface state (the upper thick line in Fig. [2j) whose wave 
function has a pronounced maximum in the region of the 
metal-organic interface (see Fig. [3]) and a surface-state- 
like penetration into the silver film. Actually, as is seen 
from Table ini), more than 50% of the IS charge accumu- 
lates in the metal film. The rest of the charge are divided 
among the interface and molecular-monolayer regions as 
~ 2/3 and 1/3, respectively. As to the energy Ejs 
of the interface state (IS) at the f point, we note that 
Eis^ which is about 0.5 eV, is quite stable with respect 
to the film thickness (see Table [I]). The dispersion of the 
interface state is adequately described by a quadratic de- 
pendence on the wave vector k with an effective mass 
mj^, which is close to that for the aforemention surface 
state for Ni > 6. 

To all appearances, the presented analysis corroborates 
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FIG. 3: (Color online) Unity- normalized charge- density distribution p/5(z,k) for the interface states as a function of z for 
different k in the ^ V ^ Xx direction, (a) The PTCDA monolayer on the nine-layer Ag(lll) film, (b) The PTCDA 
monolayer on the twelve-layer Ag(lll) film. The distributions for different /c-points shown in Fig. [2] are presented in order 
from left to right in correspondence with /c-points listed in Table |I] in order from top to bottom. 



TABLE II: Charges of the interface state in regions of the sil- 
ver film (Qb), the interface (Qi), and the PTCDA monolayer 
film (Qm) at different values of k along the X^ ^ F ^ X\ 
direction of the surface Brillouin zone (see points in Fig. [2]). 
Two cases of nine- and twelve- layer silver films are presented. 
To calculate the charges, the unity-normalized charge- density 
distribution p/5(^,k) shown in Fig. [3]has been used. 
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the interpretation of Ref. 6, in which the unoccupied dis- 
persing state of the PTCDA/ Ag(l 11) system is consid- 
ered as a genuine interface state arising from an upshift of 
the occupied Shockley surface state of the clean Ag(lll) 
surface. In fact, as our results show, there is no other un- 
occupied state, which we could associate with the state 
discussed in Refs. [6|and[5|. In our case, the only interface 
state we have found can be considered as a transforma- 
tion of the surface state that is marked as SSI in the case 
of the clean Ag(lll) films. The transformation is caused 
by the interaction of the Ag(lll) film with the PTCDA 
monolayer. 

Now, in order to establish a connection of our results 
with the STS observations of Ref. 's', we should trace the 
transformation of molecular orbitals upon moving from 
a free PTCDA molecule to the PTCDA/Ag(lll) system. 
To this end, we analyze energies, charge-density distri- 



butions p(^,k = 0), and spatial distributions of corre- 
sponding states in the molecular plane. We recognize 
electronic states in the PTCDA/ Ag(l 11) system as those 
originated from molecular levels from the localization of 
the wave function predominantly on the PTCDA mono- 
layer. Fig.|4]shows our results of such an analysis. As fol- 
lows from the figure, the LUMO level of the free molecule 
[with the structure as it appears in the PTCDA mono- 
layer on Ag(lll)] goes down under the Fermi level of the 
PTCDA/ Ag(l 11) system and splits due to the presence 
of two types of molecules (A an B) in the unit cell (see 
Fig.[T]). At that the spatial distribution, which is slightly 
modified under the influence of intermolecular interac- 
tion and adsorbate-substrate interaction, preserves prac- 
tically all its key features. Note that the wave function 
of the state with smaller energy is predominantly local- 
ized on the B-type molecule, while the state with bigger 
energy is characterized by the wave function, which is 
mainly localized on the A- type molecule. The nearly de- 
generate LUMO+1 and LUMO+2 orbitals demonstrate 
the same tendency. As well as in Refs. 15 and 19, the 
states originated from these orbitals fix their energies 
in the vicinity of leV. An analysis has shown that the 
dispersion of these states are very small as compared 
with that of the found interface state. As regards the 
spatial distribution of the latter. Fig. [5] clearly shows 
that, as well as in the case of the states arising from the 
LUMO+1 and LUMO+2 orbitals, the IS charge density 
is concentrated on the molecules. However, the charac- 
ter of the distribution on the molecule is quite differ- 
ent. Moreover, in contrast to the states arising from the 
LUMO+1 and LUMO+2 orbitals, the IS charge density 
is distributed more or less uniformly among the A- and 
B-type molecules. 

The analysis presented above does not allow us to fol- 
low the interpretation of the unoccupied dispersing state 
of the PTCDA/ Ag(l 11) system, which was made in Ref. 
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FreePTCDA PTCaVAg(lll) 

FIG. 4: (Color online) Spatial distribution of three lowest 
unoccupied orbitals of the free PTCDA molecule (at the left) 
and corresponding states in the PTCDA/Ag(lll) system at 
the r point in the case of the six- layer silver film (at the right). 




FIG. 5: (Color online) Spatial distribution of the interface 
state in the PTCDA/Ag(lll) system at the F point in the 
molecular plane for the case of the nine-layer silver film. 

[Bl, as the state formed by the molecular orbitals. At that, 
such main features of the dispersing state, which have 
been observed in the cited work by the STS, as the dom- 
inant localization on the perylene backbone and the in- 
volvement of the anhydride groups in formation of the 
state are clearly seen in Fig. [5l Here we would like to 
note that in the molecular plane the spatial distribution 
of the IS does not vary with the silver-film thickness and 
retains its form as the wave vector k changes along the 
symmetric directions of the Brillouin zone. 



IV. CONCLUSION 

In conclusion, we have investigated electronic structure 
of the PTCDA/ Ag(l 11) system within density functional 
theory with the use of the periodic slab model. The slab 
has been chosen to contain the PTCDA monolayer on the 
thin Ag(lll) film. We have considered the silver films of 
different thickness (6, 9, and 12 layers) to fix a trend in 
changes in the electronic structure of the Ag(lll) film 
upon adsorption of the PTCDA monolayer and to check 
the convergence of our results. We have found that upon 
the adsorption one of two surface states of the silver film, 
the wave function of which is localized on the clean side 
of the film, tends to have the F-point energy as that for 
the clean Ag(lll) surface. Another one transforms into 
an unoccupied interface state with the wave function that 
has a pronounced maximum in the region of the metal- 
organic interface. The interface state has the f-point 
energy of ~ 0.5 eV, which is weakly varying (within ~ 50 
meV) with the silver-film thickness. The dispersion of 
this state is nearly isotropic and adequately described by 
a quadratic dependence on the two-dimensional electron 
wave vector with the effective mass 0.43me, 0.42me and 
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OAdrUe for six-, nine-, and twelve-layer Ag(lll) film, re- 
spectively. 

We have related the aforementioned interface state 
with the unoccupied dispersing state experimentally 
observed in the PTCDA/Ag(lll) system by scanning 
tunneling^ (the effective mass of OATrrie) and two pho- 
ton photoemission^ spectroscopy (the effective mass of 
0.39 ± 0.03me). We have demonstrated that as in the 
STS observations, in the plane of the PTCDA monolayer 
film, the spatial distribution of the interface state is con- 
centrated on the molecules with quite large charge den- 
sity on the perylene backboone and anhydride groups and 
relatively small - at the perylene edges of the PTCDA 
molecule. At that, as it has been revealed by the 2PPES, 
the wave function of this state has a surface-state-like 
penetration into the silver film. We thus have corrobo- 
rated the interpretation of the origin of the unoccupied 
dispersing state of the PTCDA/ Ag(l 11) system, which 



has been proposed in Ref. 6. 

Finally, we would like to note that in our calcula- 
tion the unoccupied interface state appears at lower en- 
ergy than it is experimentally observed (0.7 eV in Ref. HI 
and 0.6 eV in Ref. 6). It can be caused by using the 
LDA to describe electronic states of the metal-organic 
interface. We believe that as in the case of the free 
PTCDA molecule (see, e.g., Ref. 20) the GW approxi- 
mation can improve the description of the electronic spec- 
trum of the system under question. 
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